A new collective phenomenon arising from spin anisotropic perturbations to a 
Heisenberg square lattice manifested in paramagnetic resonance experiments 
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We report unexpected behaviour in a family of Cu spin-| systems, in which an apparent gap in 
the low energy magneto-optical absorption spectrum opens at low temperature. This previously- 
unreported collective phenomenon arises at temperatures where the energy of the dominant exchange 
interaction exceeds the thermal energy. Simulations of the observed shifts in electron paramagnetic 
resonance spectral weight, which include spin anisotropy, reproduce this behavior yielding the mag- 
nitude of the spin anisotropy in these compounds. 

PACS numbers: 76.30-v, 75.10.Pq 



The spin (S) | two-dimensional (2D) square-lattice 
quantum Heisenberg antiferromagnet system has long 
been interesting to theoretical physicists due to the va- 
riety of transitions that can arise [U Moreover, 
the role of S = | fluctuations on a square lattice in 
the mechanism for cuprate superconductivity is hotly de- 
bated [31 IH [3]. The recently discovered family of H- 
bonded metal-organic magnets [7] offer the possibility to 
readily control the exchange parameters in a 2D system 
by changing chemical composition, thus creating spin ar- 
chitectures with desirable properties 'to order' [8] . For an 
idealized 2D system, long range magnetic order would not 
occur at finite temperature [DH]. However, in the metal- 
organic systems, interlayer coupling gives rise to a finite 
Neel temperature [Tni [TT], [T2] ■ For these quasi-2D systems 
the ordering temperature is dominated by the weakest 
(the interlayer) exchange interaction, whereas the satu- 
ration magnetic field is dominated by the strongest ex- 
change interactions, thus providing a means of estimat- 
ing the spatial exchange anisotropy in the system |5]. It 
should be noted that the more 2D the system, the wider 
the temperature (T) range, Tn < T < J/Icb, over which 
magnetic fluctuations dominate. Here we demonstrate 
that a spin anisotropy perturbation to the Heisenberg 
square lattice results in a new collective phenomenon 
within this regime, manifested as a shift in electron para- 
magnetic resonance (EPR) frequency at low temperature. 

Fig [l] illustrates the dramatic shift of the EPR mag- 
netic field and line width with temperature for a single 
frequency, /. Multi-frequency EPR measurements, fit- 
ted to / = fo + gfiBB/h, where B is the resonance field, 
show that the shifts correspond to changes in both the 
low-field intercept fo and effective g-factor g. At high T, 
fo = 0, but as the phenomenon develops with decreasing 
T, fo becomes finite and positive, behavior that strongly 
resembles an energy gap opening. Not only is this un- 
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FIG. 1: a) The Q-factor of the 26 GHz cavity-mode as a func- 
tion of magnetic field, illustrating the evolution of the EPR 
line width and resonant magnetic field with temperature in 
[CuHF2(pyz)2]C104. The data is for magnetic field applied 
perpendicular to the planes, (b) Low-temperature magneti- 
sation data in the same compound [8]. 



expected in an S — \ system; it is ruled out by the low 
temperature magnetization, which shows a monotonic in- 
crease between H = Q and the onset of saturation (see 
Fig. [TJ. Any gap would be manifested as a region of 
reduced or zero dM/dH [5]. 

We shall show below that this dramatic change is not 
linked to the antiferromagnetic transition temperature 
in either of the compounds studied, but is rather related 
to the intralayer exchange energy J. The metal-organic 
compounds [CuHF2(pyz)2]X were therefore selected for 
this study as the choice of anion molecule X can al- 
ter J controllably by a factor ~ 2 and because the J 
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values have been determined to a high accuracy using 
magnetometry [8]. In these compounds, Cu 2+ ions are 
arranged in square-lattice layers, separated by pyrazine 
(pyz) molecules. The layers are held apart by bifluoridc 
bridges [3 IE1II5]. Single crystals of [CuHF 2 (pyz) 2 ]X with 
X = CIO4 or PF 6 , were produced by an aqueous chemical 
reaction between the appropriate CuX 2 salt and stoichio- 
metric amounts of the ligands (see Q3J [14] for prepa- 
ration method details and X-ray data). Both materials 
undergo antiferromagnetic ordering, with Tn = 1.94 K 
for X = CIO4 and T N = 4.31 K for X = PF 6 . Note that 
owing to differing summing conventions the definition of 
J used in the current paper is a factor two smaller than 
that in Ref. [8], i.e. herein J p i a ne(X = CIO4) = 3.6 K 
and J planc (X = PF 6 ) = 6.2 K. 

EPR spectra of single-crystal samples were obtained in 
two ways. In the first method, the sample was placed in 
a cylindrical resonant cavity and the Q-factor and reso- 
nant frequency of the cavity were measured at each field 
point. The change in Q-factor is proportional to the mi- 
crowave absorption of the sample |15j . These measure- 
ments were carried out in the frequency range 11-40 GHz 
using a Hewlett-Packard 8722ET network analyzer. The 
cavity was placed in a 4 He flow cryostat that was capa- 
ble of stabilizing temperatures down to 1.5 K. For higher 
frequencies, the sample was placed in a confocal etalon 
resonator. The transmission, and hence change in mi- 
crowave absorption, was measured in the frequency range 
60-120 GHz using an ABmm [16 Millimetre wave Vector 
Network Analyser. The etalon was placed in a 3 He sys- 
tem that was capable of stabilizing temperatures down to 
0.6 K. For fixed temperatures, the EPR field was found 
to be a linear function of frequency, allowing the inter- 
cept and effective g-factor to be determined as described 
above; experimental values will be compared with theo- 
retical predictions below. 

Before demonstrating that the shift in intercept fo and 
effective g-factor g (Fig. 1) is a collective phenomenon 
involving a large number of spins, we show that sim- 
pler models such as dimerization cannot account for 
the data. Fig. [2^,) illustrates the Zeeman splitting for 
non-interacting electrons and b) the corresponding lin- 
ear frequency-magnetic field relationship with a zero fre- 
quency intercept [15]. Fig. [2J;) and d) includes the effect 
of an antiferromagnetic scalar exchange interaction J. 
The I TT), 73d TD + I IT)) and | ||) states will be degen- 
erate at B = 0, with the ^(| ]{) — \ {])) state separated 
by an energy — J [15] . Although this singlet-triplet gap 
causes the optically active (SS Z = il^Srot = 0) EPR 
transitions within the triplet states to become 'frozen out' 
at T < J, it does not introduce a gap across which opti- 
cally active EPR transitions occur, which would lead to 
a finite frequency intercept. 

Fig. [2^) and f ) includes a spin anisotropy in the ex- 
change interaction, J zz ^ J xx = J yy , lifting the zero field 
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FIG. 2: Illustration of how a spin-anisotropic exchange in- 
teraction can give rise to a an EPR spectrum with a finite 
frequency intercept in the case of an antiferromagnetically 
coupled dimer. Respectively (a), (c) and (e) show the Zee- 
man splitting for a non-interacting, isotropically coupled and 
spin-anisotropically coupled dimer. (b), (d) and (f) illustrate 
the corresponding frequency-magnetic field relationships of 
the resonant absorption. 



degeneracy of the ST 0t = 1, S x = ±1 and S^ot = l,S z = 
states (Fig. [2]), which in turn lifts the field degeneracy 
of the two intratriplet transitions, separating EPR lines 
2 and 3. Although some of the possible EPR lines now 
have finite frequency B — intercepts, this simple model 
does not reproduce the systematics of the effect reported 
here. ^(| ][) — \ [])) is the groundstate so that the 
two transitions potentially observable at low temperature 
will be lines 1 and 4; however, for these to be optically- 
active, off-diagonal exchange terms (e.g. J xy ) such as the 
Dzyaloshinsky-Moriya interaction must be present (since 
they mix the singlet and triplet states) |15j . We note 
that despite a Dzyaloshinsky-Moriya term arising from 
a lower order spin-orbit perturbation than the diagonal 
spin anisotropy terms, that they are precluded by the in- 
version symmetry about the mid point of the dominant 
exchange interaction. At elevated temperatures T > J, 
one would observe an EPR spectrum dominated by lines 
2 and 3 (each with a finite B = intercept) ; as the sam- 
ple cooled to T < J, this would change to a spectrum 
dominated by the much weaker (or even absent) lines 1 
and 4, with an interecept differing from that of 2 and 3. 
Such behavior does not lead to a smooth thermal evolu- 
tion of the EPR line or the preservation of spectral weight 
to temperatures T < J that we observe. 

Having shown that a simple local distortion cannot 
cause the effect shown in Fig. 1, we now turn to a finite 
cluster approach applied to an anisotropic 2D Heisenberg 
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model with spin anisotropy for two field orientations: 
H = JJ2 S'Sj" + A(SfSf + SfSf) - g Tt i B B £ S? < y , 

where A is the spin anisotropy, with A = 1.0 being the 
isotropic (pure Heisenberg) case, A < 1 an Ising-like or 
easy-axis antiferromagnet and A > 1 an easy-plane xy- 
antiferromagnet. It should be noted that the input pa- 
rameter to the model gr is a temperature independent g- 
tensor that reproduces the high temperature (T >> J) g- 
factor anisotropy arising from the spin-orbit interaction. 
This g- anisotropy is consistent with the magnetic d x 2_ y 2 
orbital lying in the 2D planes [5] . J pcr p is assumed to be 
zero, since J p0 rp << -/plane; a s demonstrated in [S . This 
model was calculated using full diagonalization at finite 
temperature for 12, 14 and 16 sites [T7]. There was little 
difference between the results with different numbers of 
sites, and therefore only the 14 site data are displayed 
here. The simulations were carried out for anisotropy 
values A=1.02, 1.04, 1.06 and 1.08. It was found that in- 
creased spin anisotropy leads to the parallel shift of the 
peak away from its initial position while at the same time 
the peak gains a finite width. Multiple smaller peaks in 
addition to the main peak were observed in the simula- 
tions for A ^ 1, due to finite-size effects. 

To compare the simulations to the results for 
[CuHF 2 (pyz) 2 ]C104, the resonant magnetic field, Lorentz 
linewidth and spectral weight of the resonance were cal- 
culated for a resonance at 26 GHz. The resonant mag- 
netic field was calculated as -Br, = 2Bq — (u), where 
B = 0.35J (26 GHz in units of J) and 

_ / LL>S xx (u),B )du> 
J S xx (u),B )dui 

The Lorentz linewidth 

V Kel) 

where a-> rc i = lu — (to) . The spectral weight was calculated 
as: 

/= J VVo,5)d£~(l-e-^ T ) Js xx (tu,B )du 

with the dynamic spin structure factor 
1 f°° 

S xx (lu) = -Re / e^S^S^O)^. 

A Lorentzian fit to the experimental data yielded the 
linewidth and spectral weight. As can be seen from 
Fig. [3] the results of the simulations reproduce the salient 
features of experimental results for [CuHF2(pyz)2]C104. 
The experimental values clearly lie between the simu- 
lated values of A = 1.04 and 1.08. Taking into account 
the fact that above T = 0.8 J the 1.04 simulation pro- 
vides a good fit for both the position of the resonance 



and its linewidth, the best match for the data is given by 
a spin anisotropy of 1.04. 

From Fig [3] the changes in the position and width 
of the EPR line in [CuHF 2 (pyz) 2 ]ClC>4, start as the 
temperature is lowered through 3.5 K; this corresponds 
to T/J=l (as opposed to T N = 1.94 K). Turning to 
[CuHF 2 (pyz) 2 ]PF 6 , Fig. [2] shows that both g and / 
undergo a dramatic change on cooling through 4.5 K. 
Coincidcntally, this is close to the ordering temperature 
Tn [H] . The greater relative thermal separation of Tn & n d 
T = J in the CIO4 compound than in the PFg salt re- 
flects the smaller interlayer exchange energy; one might 
say that the former is a closer approximation to two di- 
mensionality than the PFq compound |8J . In spite of the 
proximity of the ordering temperature to the onset of 
the EPR shifts, the fit of the model to the data from the 
PFg compound that we will now give shows that it is the 
exchange interaction J that determines the temperature 
scale of the effect and not Tn- 

To compare the simulations for the results for 
[CuHF 2 (pyz) 2 ]PF 6 the experimental g-factors and fre- 
quency intercepts were calculated. The theoretical g- 
factors and intercepts are given by the slopes and in- 
tercepts of the ui, B) plots, which were obtained from 
the extrapolation through five calculated points, uj/3 — 
1.0,1.1,1.2,1.3,1.4 for correspondence with the experi- 
mental frequency range. 

In Fig. [4] we can see that the behaviour of the exper- 
imental g-factor and intercept are well matched by the 
simulations, with the dramatic difference between field B 
perpendicular to and B parallel to the square Cu plane 
being reproduced. For the B parallel data, the fact that 
we do not observe an upturn in the g-factor and the small 
positive intercept at low temperature is most likely due to 
a small angular misalignment (since the B perpendicular 
effects are so much larger). 

Considering the B perpendicular data, the g-factor ex- 
perimental data suggests a spin anisotropy above 1.08, 
whereas the intercept data suggests an anisotropy close 
to 1.04. If we consider the corresponding data for X 
= CIO4 (for which data was taken at only a few values 
due to the extreme weakness of the signal at low tem- 
peratures) we find that at T/J — 0.27 the experimental 
values are g = 0.96 go and intercept = 0.12 J. This gives 
a similar anisotropy, 1.06 - 1.04. For both materials it 
can be clearly seen that for T/J < 1 the experimental 
changes in both intercept and g factor are substantially 
more rapid than the simulations suggest. We propose 
that this is due to the experimental temperature scale 
being contracted relative to the theoretical temperature 
scale upon entering the ordered state, i.e. the phase tran- 
sition has the effect of contracting the entropic (theoreti- 
cal temperature) scale. The shape of the intercept curve 
suggests that the experimental value of the intercept has 
nearly stabilised by the lowest T/J values and therefore 
the intercept is taken to give the best indication of the 
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FIG. 3: A comparison of the experimental EPR data, cir- 
cles, from [CuHF2(pyz)2]C104 measured at 26 GHz with the- 
oretical simulations, A =1.02, 1.04, 1.06, 1.08, 1.10; solid, 
dot, dash, dot-dash and dot-dot-dash respectively. The left 
and bottom axes are absolute units, the right and top axes 
are renormalized by the exchange energy and spectroscopic 
g-factor. a) The variation of the magnetic field at which reso- 
nance occurs as a function of temperature, b) The linewidth 
(full width at half maximum of a Lorentzian fit) of the res- 
onance peak as a function of temperature, c) The spectral 
weight (area under Lorentzian fit) normalized to its 15 K 
value. In all the cases error bars are smaller than data points. 
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FIG. 4: Comparison of the experimental g-factor for X = 
PFa to simulations for different values of the spin anisotropy, 
A =1.02, 1.04, 1.06, and 1.08, solid, dotted, dashed and dot- 
dashed lines respectively, for (a) B perpendicular to and (b) 
B parallel to the Cu square latttice. Also, comparison of 
the experimental frequency intercept of resonant frequency 
vs magnetic field and simulated values for (c) B parallel and 
(d) B perpendicular. 



environment in which to characterize this effect, strongly 
coupled copper octahedra are ubiquitous in correlated 
electron systems; for example, J ~ 100's K in the parent 
phase of the high Tq superconductors. 
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spin anisotropy (around 1.04 in both compounds). The 
approximate conservation of spectral weight in both ex- 
perimental and theoretical results also demonstrates the 
robustness of our model. 

In conclusion we have observed a collective phe- 
nomenon in two members of the organic magnet system 
(Cu(HF2)(pyz)2)X that produces a shift in the frequency 
intercept of the EPR data which resembles, but does 
not correspond to, a gap opening in the system. Analo- 
gously to an anisotropic g-tensor, the spin anisotropy in 
the exchange interaction responsible for this effect most 
likely originates from spin-orbit coupling. As a result, 
this collective phenomenon is expected to be strongest 
for low-spin transition metal ions in relatively low sym- 
metry environments, like the octahedral copper site (3d 9 ) 
in (Cu(HF2)(pyz)2)X. Although the organic systems in- 
vestigated provide an ideal 'low exchange energy scale' 



* Electronic address: rmcd@lanl.gov 

[1] M. Kastner, R. Birgenau, and G. S. amd Y. Endoh, Rev. 
Mod. Phys. 70, 897 (1998). 

[2] U. Schollwock, D. Farnell, and R. Bishop, eds., Quantum 
magnetism (Springer, Berlin, 2004). 

[3] E. Manousakis, Rev. Mod. Phys. 53, 1 (1991). 

[4] J. Schrieffer and J. Brooks, eds., High temperature su- 
perconductivity theory and experiment (Springer, Berlin, 
2007). 

[5] S. Julian and M. Norman, Nature 447, 537 (2007). 

[6] R. McDonald, N.Harrison, and J. Singleton, J. Phys. CM 
21, 012201 (2009). 

[7] J.L. Manson, M.M. Conner, J.A. Schlueter, T. Lan- 
caster, S.J. Blundell, M.L. Brooks, T. Papageorgiou, 
A.D. Bianchi, J. Wosnitza, M.H. Wangbo, Chem. Com- 
mum. p. 4894 (2006). 

[8] P.A. Goddard, J. Singleton, P. Sengupta, R.D. McDon- 
ald, T. Lancaster, S.J. Blundell, F.L. Pratt, S. Cox, N. 
Harrison, J.L. Manson, H.I. Southerland, J.A. Schlueter, 
New Journal of Physics 10, 083025 (2008). 



[9] N. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 
(1966). 

[10] N. Christensen, H. Ronnow, D. McMorrow, A. Harri- 
son, T. Perring, M. Enderle, R. Coldea, L. Regnault, 
and G. Aeppli, Proc. Nat. Acad. Sci. 104, 15264 (2007). 

[11] J. Choi, J. Woodward, J. Musfeldt, C. Landee, and 
M. Turnbull, Chem. of Materials 15, 2797 (2003). 

[12] T. Lancaster, S. Blundell, M. Brooks, P. Baker, F. Pratt, 
J. Manson, M. Conner, F. Xiao, C. Landee, F. Chaves, 
et al., Phys. Rev. B 75, 094421 (2007). 

[13] J.L. Manson, H. Southerland, J. Schlueter and K. Funk, 
preprint (2008). 

[14] T. Lancaster, S.J. Blundell, P.J. Baker, M.L. Brooks, 
W. Hayes, F.L. Pratt, J.L. Manson M.M. Conner, J. A. 
Schlueter, Phys. Rev. Lett. 99, 267601 (2007). 

[15] A. Abragam and B. Bleaney, Electron Paramagnetic Res- 
onance of Transition Ions (Oxford University Press, New 
York, 1970). 

[16] R. D. McDonald et al., Rev. Sci. Instrum. 77, 084702 
(2006). 

[17] for B//y the maximum number of sites used was 14 due 
to the lack of S z symmetry. 



